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Interference cancellation method and receiver 

field of the invention 

[0001] The invention relates to a receiver and a method of 
suppressing a narrow-band interference from a wide-band signal. 

BACKGROUND OF THE INVENTION 

[0002] Telecommunication connections are subject to various 
interferences. There are several reasons for interference. In 
telecommunication connections, the transmission path used for transmitting 
signals is known to cause interference to telecommunication. This occurs 
regardless of the physical form of the transmission path; whether the 
transmission path is, for instance, a radio link, an optical fibre or a copper 
cable. Particularly in radio telecommunication there are frequently situations 
where the quality of the transmission path varies from one connection to 
another and also during a connection. 

[0003] Radio path fading is a typical phenomenon that causes 
changes in a transmission channel. Other simultaneous connections may also 
cause interferences and they can vary as a function of time and place. 

[0004] To reduce effects of various interferences, several solutions 
have been evolved. This is the case in radio communication in particular. 
Various telecommunication methods have been evolved, aiming at achieving 
as high spectral efficiency as possible and still minimising the effect of 
interference. A wireless telecommunication method which has been very much 
studied lately is a wide-band spread-spectrum multiple access method CDMA 
(Code Division Multiple Access). 

[0005] Compared to multiple access methods that were generally 
used earlier, the CDMA has a plurality of advantages. In a CDMA method, a 
narrow-band data signal of the user is multiplied by a spreading code having a 
considerably wider band to a relatively wide band. Bandwidths used in known 
experimental systems include 1.25 MHz, 10 MHz and 25 MHz. In the 
multiplying process, the data signal spreads to the entire band used. All users 
transmit simultaneously by using the same frequency band. Each connection 
between a base station and a mobile station uses its own spreading code, and 
the signals of the users can be distinguished from each other in the receivers 
on the basis of the spreading code of each user. The aim is to select the 
spreading codes so that they are mutually orthogonal, i.e. they do not correlate 
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with each other. The above described CDMA method is called a direct- 
sequence method, DS - CDMA (Direct Sequence CDMA). There are other 
CDMA methods as well, such as a frequency hopping method FH-CDMA 
(Frequency Hopping CDMA), in which the frequency used is varied quickly 
5 according to the used spreading code. A combination of these methods is also 
possible. 

[0006] A problem of wide-band data transmission is narrow-band 
interference, which is typically caused by narrow-band signal sources external 
to the system, using the same or overlaying frequency band with the system. 

10 Typical of these interference signals is that their properties and structure often 
differ from the signals of the system considerably. They often use a lower data 
transmission rate. Also, it is often so that the properties of interfering signals 
are not known in advance. 

[0007] The CDMA can, to some extent, inherently compensate the 

1 5 effect of narrow-band interference on wide-band data transmission, but if the 
interfering signal is much stronger than the wide-band signal, it may cause 
considerable interference. 

[0008] To suppress narrow-band interference from a wide-band 
signal, various methods have been evolved. The publication S. Sandberg, 

20 "Adapted Demodulation for Spread Spectrum Receivers which Employ 

Transform-Domain Interference Rejection", IEEE Trans. On Communications, 
Vol. 43, pp 2502-2510, Sept. 1995, discloses a method in which the problem is 
approached in the frequency domain, and the interference is suppressed by 
removing the frequencies where the interference is assumed to be. In the 

25 presented solution, however, the location of the interference in the frequency 
band and the bandwidth of the interference are assumed to be known. The 
interference is also assumed to be stationary. 

[0009] The publication M. Lops, A. Tulino, "Automatic Suppression 
of Narrowband Interference in Direct Sequency Spread-Spectrum Systems", 

30 IEEE Trans. On Communications, Vol. 47, pp. 1133- 1136, August 1999, 
discloses a method in which non-stationary interference, the properties of 
which are not known, can be suppressed. However, the achieved results are 
not qualitatively satisfactory. 
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BRIEF DESCRIPTION OF THE INVENTION 

[0010] It is an object of the invention to implenaent an improved 
method and equipment implementing the method for interference cancellation. 
This is achieved by a method of suppressing narrow-band interference from a 
5 wide-band signal, wherein a signal is received, an orthonormal conversion of 
the signal into subspace components of a desired subspace Is performed, the 
converted signal is decoded with a decoder, whereby an estimate for the 
received signal is obtained, and when an estimate for narrow-band 
interference properties is determined, the estimate obtained from an output of 

10 the decoder is subtracted from the received signal before the orthonormal 
conversion is performed, and, by means of the determined estimate, effect of 
the subspace components comprising narrow-band interference signals is 
reduced in the received signal. 

[0011] The invention also relates to a receiver comprising 

15 interference cancellation means for suppressing a narrow-band interference 
signal from a received signal, means for performing an orthonormal 
conversion of the signal into subspace components of a desired subspace, a 
decoder connected operationally to the output of the interference suppression 
means, in which decoder an estimate for the received signal is obtained. In the 

20 receiver of the invention, the output of the decoder is operationally connected 
to the interference suppression means, and the conversion means are 
arranged to determine an estimate for narrow-band interference properties, in 
which determination the estimate obtained from the output of the decoder is 
subtracted from the received signal before the orthonormal conversion is 

25 performed, and by using the determined estimate, the interference 
cancellation means are arranged to reduce effect of the subspace 
components comprising narrow-band interference signals in the received 
signal. 

[0012] The preferred embodiments of the invention are disclosed in 
30 the dependent claims. 

[0013] The invention is based on projecting the received signal onto 
an orthonormal subspace. A desired signal has a projection deviating from 
zero onto most orthonormal bases, whereas narrow-band interference has 
concentrated on specific directions. By suitably selecting the orthonormal base 
35 and the projection to be performed, the interference can thus be reduced or 
suppressed in its entirety. By subtracting the preliminary estimate of the 
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desired signal from tiie received signal, the subspace components on which 
the interference concentrates can be estimated better than before. Thus, the 
projection onto the orthonormal components orthogonal to the Interference 
can be performed better than before. The method can be repeated iteratively 
5 many times. 

BRIEF DESCRIPTION OF THE FIGURES 

[0014] In the following the invention will be described in greater 
detail in connection with the preferred embodiments, with reference to the 
attached drawings, in which 
10 [0015] Figure 1 shows an example of a system according to an 

embodiment of the invention, 

[0016] Figure 2 shows a second example of a system according to 
an embodiment of the invention, and 

[0017] Figure 3 illustrates a system model according to an 
1 5 embodiment of the invention, 

[0018] Figure 4 illustrates a structure of a receiver according to an 
embodiment of the invention, and 

[0019] Figures 5a and 5b illustrate solutions according to preferred 
embodiments of the invention. 

20 DESCRIPTION OF THE PREFERRED EMBODIMENTS OF THE INVENTION 

[0020] Preferred embodiments of the invention can be applied in 
telecommunication systems where a desired signal is transmitted and received 
as a wide-band one. Such a telecommunication system is a wide-band CDMA 
radio system, for instance. In the following example, preferred embodiments of 

25 the invention will be described in a universal mobile telephone system 

employing a wide-band code-division multiple access method, yet without 
restricting the invention thereto, 

[0021] With reference to Figure 1 , a structure of a mobile telephone 
system is explained by way of example. The main parts of the mobile 

30 telephone system are core network CN, terrestrial radio access network of the 
mobile telephone system UTRAN (UMTS terrestrial radio access network) and 
user equipment UE. The Interface between the CN and the UTRAN is called lu 
and the air interface between the UTRAN and the UE is called Uu. 

[0022] The UTRAN comprises radio network subsystems RNS. The 

35 interface between the RNSs is called lur. The RNS comprises a radio network 
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controller RNC and one or more nodes B. The interface between the RNC and 
B is called lub. The coverage area, or cell, of the node B is marked with C in 
the figure. 

[0023] The description of Figure 1 is relatively general, and so it is 
5 clarified with a more specific example of a cellular radio system shown in 

Figure 2. Figure 2 only includes the most essential blocks, but it is obvious to a 
person skilled in the art that the conventional cellular radio network also 
includes other functions and structures, which need not be further explained 
herein. It is also to be noted that Figure 2 only shows one exemplified 
10 structure. In systems according to the invention, details can be different from 
what are shown in Figure 2, but as to the invention, these differences are not 
relevant. 

[0024] A cellular radio network thus typically comprises a fixed 
network infrastructure, i.e. a network part 200, and user equipment 202, which 

15 may be fixedly located, vehicle-mounted or portable terminals. The network 
part 200 comprises base stations 204. A base station corresponds to the node 
B shown in the previous figure. A plural number of base stations 204 are, in 
turn, controlled in a centralised manner by a radio network controller 206 
communicating with them. The base station 204 comprises transceivers 408 

20 and a multiplexer unit 212. 

[0025] The base station 204 further comprises a control unit 210 
which controls the operation of the transceivers 208 and the multiplexer 212. 
The multiplexer 212 arranges the traffic and control channels used by several 
transceivers 208 to a single transmission connection 214. The transmission 

25 connection 214 forms an interface lub. 

[0026] The transceivers 208 of the base station 204 are connected 
to an antenna unit 218 which is used for implementing a bi-directional radio 
connection 216 to the user equipment 202. The structure of the frames to be 
transmitted in the bi-directional radio connection 216 is defined separately in 

30 each system, the connection being referred to as an air interface Uu. 

[0027] The radio network controller 206 comprises a group 
switching field 220 and a control unit 222. The group switching field 220 is 
used for connecting speech and data and for combining signalling circuits. The 
base station 204 and the radio network controller 206 form a radio network 

35 subsystem 224 which further comprises a transcoder 226. The transcoder 226 
is usually located as close to a mobile services switching centre 228 as 
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possible, because speech can then be transferred in a cellular radio network 
form between the transcoder 226 and the radio network controller 206, which 
saves transmission capacity. 

[0028] The transcoder 226 converts different digital speech coding 
5 forms used between a public switched telephone network and a radio network 
to make them compatible, for instance from a fixed network form to another 
cellular radio network form, and vice versa. The control unit 222 performs call 
control, mobility management, collection of statistical data and signalling. 

[0029] Figure 2 further shows the mobile services switching centre 
10 228 and a gateway mobile services switching centre 230 which controls the 
connections from the mobile communications system to the outside world, in 
this case to a public switched telephone network 232. 

[0030] The solution according to the preferred embodiments of the 
invention can be applied to both a base station receiver and a user equipment 
1 5 receiver. 

[0031] Let us examine the example of a system model shown in 
Figure 3. The model is somewhat simplified, and it does not include, for 
instance, radio frequency parts, antennas, chip pulse form scramblers and A/D 
converters at transmission and reception ends, which are typically used in 

20 radio systems. It is assumed in this example that the system uses both a 
frequency hopping (FH) technique and a direct-sequence (DS) technique. 
However, it is to be noted that the embodiments of the invention are not 
restricted to this combination. Let us assume that binary data flow d(k) 300 to 
be transmitted is coded by a turbo-coder 302 with a code rate 1/R by using Ne 

25 parallel recursive systematic coders, each of which uses a code rate 1/Re and 
which are separated from each other by interleavers. Code bits are mapped 
with a BPSK method, and the symbol length is Tb. Each data symbol b(kR+j), j 
= 0,...,R-1 is then modulated in an DS-BPSK modulator 304 with a spreading 
code p(t) 306 which is generated in a spreading code generator 308. The 

30 complex envelope of a useful signal 310 before a FH modulator 310 is 

ST(t) = Z S b(kR + j) p(t-kRTt>-C-1)Tb) 

= Z Z b( kR + j ) X QUtc( t - iTc - kRNTc - 0-1 )NTc ), 

^=0 J=Q i=0 
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where Q is the number of infornnation bits per one frame, [b(kR), b(kR + 1),..., 
b(kR + R-1)] = Bk, where k = 0,...,Q-1 is the code word which is transmitted in 
the k;th information signalling time slot RTb, Tc is the chip time slot (i.e, the bit 
time slot of the spreading code), and which is proportioned to the signalling 
5 time slot Tb such that Tb = NTc, where N is processing gain, c = [Co,...,Cn-i]^ is 
the spreading sequence and u(t) is defined as follows: 

u(.)4^' 

0, elsewhere 

[0032] A DS signal 310 Is modulated in the FH modulator 312 
where frequency hops are generated according to a pseudo-random code 316 
10 generated in the FH generator 314. 

[0033] The signal is transmitted to a channel 318 where 
interference and noise 320 are added to the signal. 

[0034] The signal that has passed through the channel is received 
in the receiver, amplified in the radio frequency parts and passed to the FH 
1 5 demodulator 322 where frequency hopping is decoded according to the 
pseudo-random code 326 generated in the FH generator 324. 

[0035] Let us next examine the received signal in the output of the 
FH demodulator 322. Then the signal has the form 

r(t) = Aej%(t-T) + i(t) + w(t), (1) 
20 where Ae^"^ is the complex gain caused by the channel effect, x is the random 
delay caused by the channel and w(t) is thermal noise which can be modelled 
as a complex Gaussian process with a power spectral density 2No. i(t), for its 
part, is the complex envelope of the narrow-band interference signal. 

[0036] Let us assume that the delay t of the desired signal is either 
25 known or it is estimated. Then, the desired signal can be expressed in the 1 = 
(k+j):th signalling time slot of the code bit as a vector in an N-dimensional 
subspace from those functions in the space L^ [1 Tb + t, (l +1 )Tb + t[ whose 
quadratic absolute value can be integrated over the duration of Tb. In the 
following, the marking Sn(l) and the term signal space are used for this 
30 subspace. 

[0037] When the received waveform (1) is projected onto the 
selected orthonormal base and the entire frame is observed, the following 
NRQ-dimensional observation vector is obtained: 
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r = [r^(1),...,r''(l),...T'"(NRQ)f, (2) 

r(l) = Ae^%(l)s + i(l) + w(l). 
s, i ja w indicate N-dimensional vectors of the projections of the desired signal, 
interference and noise onto a suitable orthonormal base b of the signal space 
5 Sn(l). s is proportional to spreading codes c. In the above, refers to 
transpose. 

[0038] An idea of the invention is that since the desired signal is of 
a spread-spectrum nature and the noise is uncorrelated, it is possible to find 
an orthonormal conversion, whereby the desired signal has spread evenly to 

10 all N components, whereas the narrow-band interference has concentrated on 
a subspace of the signal space, having small dimension. When a strong 
narrow-band interference has been received, some of the terms of r have an 
absolute value which is considerably higher than that of others and by which 
they can be separated from the others, and thus interference cancellation can 

15 be formed as a problem of finding these unknown components. When these 
components have been found, their effect can be reduced either by excising 
them to the desired level or by nullifying components. One embodiment is to 
include only those samples in further processing whose absolute value is 
below the given limit. This alternative is called type 1 censoring. Another 

20 embodiment is to nullify L samples having the highest absolute values, where 
L is the estimate of the subspace dimension. This alternative is called type II 
censoring. 

[0039] From the output of the FH demodulator 322, the signal is 
thus passed to interference cancellation means 328 which perform the above 

25 described projection. After interference cancellation, DS-BPSK demodulation, 
i.e. the decomposition decoding 330 can be performed with a spreading code 
p(t) 332, which is generated in a spreading code generator 334. The decoded 
signal is passed to the turbo-decoder 336. 

[0040] Another idea of the invention is that since the desired signal 

30 may also have components in the same direction as the narrow-band 

interference, a more accurate estimate for the interference is obtained, when 
the estimate of the desired signal is first subtracted from the received signal, 
whereby the signal to be processed is such that the portion of interference and 
noise can be detected more easily. The subspace components of the 

35 interference in this residual signal can be estimated more accurately. Thus, 
the signal is fed back from the turbo-decoder 336 to the interference 
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cancellation 328. The estimated signal is multiplied by the spreading code in a 
multiplier 338, before it is taken to the interference cancellation. 

[0041] Let us examine an example of the implementation according 
to a preferred embodiment of the invention by means of Figure 4. The figure 
shows two blocks, an interference cancellation block 400 and a decoder block 
402. In the preferred embodiment of the invention, the decoder block is 
implemented as a turbo-decoder. A decoder block of another type is possible, 
too. From the previous iteration round, the received signal r(l) 404 and the 
feedback soft decisions b^ 406 are supplied as input to the interference 
cancellation block from the output of the turbo-decoder 402. These soft 
decisions are subtracted from the received signal. Before the subtraction, they 
must first be multiplied by the spreading code into a direct-sequence spread- 
spectnjm form. On each Iteration round, the interference cancellation block 
first calculates an orthonormai conversion U'', estimates the Interference 
subspace SfCi) (i.e. the main direction for S|(l)), and finally projects the 
signal onto the subspace 8^(1) orthogonal to the interference subspace. 

[0042] Let us examine the implementation of an orthonormai 
conversion a bit later. Let the term z^(l) refer to the absolute value vector of 
the vector f = [r ""(I),..., f^{i),..., r''(NRQ)]"^, where 

r(l) = U^(r(l)-Ae'*bP(l)s) (3) 

and let z~(i) refer to its ranked version, i.e. 

Ai) = [zUl) z^_i(i)f = [| fo^^d)! \f^, (i)|f 

zE(i) = [z[;)(l) Z^N) (Df, 

where z^^ls the i:th grade statistics of 2^(1), and b^(l) = sign [ m { L^(l)}] is 
the hard estimate of the code bit when L^ (1) Is the soft output decoded from 
the interleaving of the turbo-decoder. According to formula (3), the 
estimates'^ (1) of the desired signal, multiplied by the spreading code s, is thus 
subtracted from the received signal r(l). 

[0043] When L samples having the highest absolute values are 
nullified, we have a new data vector v''(i): 

v^(l) = G'^(l)U^(l), (4) 
where the censoring matrix G^{1) Is an N-dimensional quadratic diagonal 
matrix which depends on the received vector r(l) and Is defined as follows: 
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type I 



G^(l)i.i 



1, if D,(l)>z^,(l) 
0, if D,(l)<z|;)(i) 



(5) 



type II 



G"(l)i,i 



1. if VdO^TO) 
0, if Z(,.„(l)< r;^(|) 



where Di(l) indicates the suitable threshold values and L is the estimate of the 
interference subspace Sf (l) dimension. The superscript O"^ indicates that the 
value of the term Is from the previous iteration round. Type I and type II thus 
refer to the earlier mentioned censoring methods of suppressing interference. 
Naturally on the first iteration round, no estimation results of the previous 
iteration round are available, and then b^ = 0 can be used as initial values. 

[0044] If the calculation needs to be made less complex, instead of 
soft decisions obtained from the decoder, hard decisions given by the decoder 
can also be used in calculating the censoring matrix G^(i). Then the censoring 
matrix is defined as follows: 



integer, which is not bigger than x. 

[0045] Let us next examine the selection of the parameters, such 
as D and L, which are to be used in the calculation. It is also essential for the 
calculation to define the base b of the orthonormal conversion. The 
interference cancellation Is most effective. If the projection of a narrow-band 
interference signal in the signal space comprises a limited number of 
components deviating from zero. Therefore, the conversion base should be 
selected so that it focuses the narrow-band interference on a subspace in the 
signal space, having small dimension. 

[0046] In a preferred embodiment of the invention, the covariance 
matrix of the vector r^{i) is estimated, where r^{i) is based on N- 

dimensional data vectors r = [r(0) r(RQ-1)], which are received during one 

frame, i.e. 




(6) 



where I — J is in accordance with formula (5) and LxJ is the biggest 
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M^P,P=^E"V>)fP(l) (7) 

[0047] By using eigenvalue value decomposition EVD, the 
orthonormal conversion can be defined easily from the formula 

5 where A denotes the singular value matrix. When this method is employed, on 
each iteration round a sample covariance matrix r ^(1) has to be calculated 
and its eigenvalue decomposition has to be evaluated. 

[0048] Another preferred embodiment of the invention utilizes the 
assumption of the interference having a narrow band, i.e. that its bandwidth is 
10 substantially smaller than the bandwidth of the desired signal. The Fourier 
base can thus be selected as a base, i.e. 

^1 (t) = UT^(t-mTJexpH^^ 1 = 0,...,N-1. (8) 

[0049] In this case, the orthonormal conversion corresponds to an 
N-point DFT conversion of the signal sampled from the output of the pulse 

15 form scrambler. When the received signal is corrupted by interference with a 
few frequency peaks only, this solution is particularly advantageous. 

[0050] In a preferred embodiment of the invention, it is also 
essential to determine the threshold values Dj and the interference space 
dimension L mentioned in formula (5). If the eigenvalues obtained from 

20 formula (7) are evaluated and a suitable information criterion, such as so- 
called Akaik (AlC), the smallest description length or an eigenvalue threshold 
method, is used, the interference subspace Sj (1) dimension can be estimated 
adaptively. These information criteria are prior art methods, and the smallest 
mapping length, for instance, is explained in the publication M.Wax, T.Kailath: 

25 "Detection of signals by information theroretic criteria", IEEE Trans. Acoust. 
Speech, Signal Processing, Vol. 33, pp. 387-392, 1985, which is incorporated 
herein as a reference. The parameter L can thus be set according to the 
dimension, 

[0051] The threshold value Di, for its part, can in principle be set so 
30 that it maximises the value 

PriMl)i>Dm|i^O). 
i.e. the detection probability, provided tliat the false alarm probability 
Pr{|r^(l)|>D, |i = 0} 
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Var 



^0 

i2 



is fixed to a desired level. This maximisation can be calculated by means of a 
known likelihood test. This method requires, however, the dependence on the 
parameters of the interference signal, and it is not optimal. In a preferred 
embodiment of the invention, a more intuitive solution model is applied, where 
only the statistics of the term \r^il)\ is utilized and a situation free of 
interference is assumed. When a threshold is set, it is essential to select them 
in such a manner that the signal is not censored unnecessarily when the 
interference is weak and yet so that the Interference cancellation is ensured 
when the interference is strong. Let us assume that the eigenvalue E[b (l) 
w(i)| i=0] = 0; as a result, when the conditions of b (i) of b(l) and when i=0 
and when |rm|^ is the non-centralised variable of the -i distribution, having two 
degrees of freedom, then 

= A^|s„|^E[(1-b(l))2] +2Nc 
= 4No +4NoAi E[(1-b(l))2]|sr 

where e[(1- 6(1 ))2] = (2-2E[b(l)b(l)]) = (^2-2sign[A^(/)]tanh |[a^(I)' 

[0052] Let us examine the example of the implementation 
according to the preferred embodiment of the invention, which was already 
partly explained and shown in Figure 4. The figure thus shows two blocks, the 
interference cancellation block 400 and the decoder block 402 which is 
preferably, but not necessarily, a turbo-decoder. From the previous iteration 
round, the received signal 404 and the feedback soft decisions 406 are 
supplied as input to the interference cancellation block 400 from the output of 
the turbo-decoder 402. On each iteration round, an upper calculating unit 408 
of the interference cancellation block 400 first determines an estimate for the 
narrow-band interference properties by subtracting the estimate obtained from 
the output of the turbo-decoder from the received signal. Then an orthonormal 
conversion is calculated on the basis of formulas (7) or (8), for instance, 
and the interference subspace Sf (1) is estimated from the residual signal. 
Then the upper calculating unit defines the censoring matrix G'^on the basis of 
formula (5), for instance. 

[0053] The information on the orthonormal conversion and 
censoring matrix is transmitted to a lower calculating unit 410 where the actual 
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censoring is performed and tlie new data vector v''(i) according to formula (4) 
is calculated. This data vector is passed to the turbo-decoder 402 as input. 

[0054] Both calculating means 408, 410 of the interference 
cancellation block 400 can preferably be implemented programmatically by 
5 means of a processor and a suitable program or also as an ASIC circuit or by 
means of separate logic circuits. It is to be noted that the interference 
cancellation block can also be implemented by means of one or more 
calculating units and the implementation described herein is only an example, 
as is obvious to a person skilled in the art. 

10 [0055] The turbo-decoder 402 comprises a SISO (soft-in soft-out) 

decoder. The turbo-decoder comprises two parallel-concatenated decoders 
412, 414, between which interleaving 416, 418 is performed. The decoding is 
based on decoding the component codes alternately and transmitting so- 
called extrinsic information, which is a part of the soft output of the SISO 

15 decoder, to the next decoding stage. The signal can be iterated several times 
in the turbo-decoder by re-feeding the obtained estimates 420 to the decoder. 
The deinterleaving 422 of is performed therebetween. The output of the turbo- 
decoder is fed back 406 to the interference cancellation block, and the loop 
also includes the deinterleaving 424. The decoder can also comprise a control 

20 block 426 controlling the function of different parts of the decoder. If several 
iteration rounds are perfomned in the decoder, the control block 426 gives a 
command to stop the iteration and transmits the estimates further to the 
interference cancellation block. The estimates are also taken further to other 
parts 428 of the receiver. The iteration can also be controlled from outside the 

25 decoder at a higher level of the receiver. 

[0056] On the P:th iteration round, the decoder makes a decision 
on the transmitted bits on the basis of the vector v'', which has the form 

v^ = [v^^^l ),..., v^^ (1) vP>RQ)]^ = G^'u'^r, 

where v^ (l) is obained from fonnula (4) and G^U^ is an NRQ-dimensional 

30 matrix which has the form 

G^U'' = diag ( G^(0)U^...,G^(RQ)U''). 
For the sake of clarity, the term 

refers in the following to an NRe-dimensional vector which has the form 
35 (k) = [vf^"'(k),v'^^(k+hRe+1 ) v''^(k+(h+1 )Re-1 )f. 
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h = 0....,Ne-1 k = 0,...,Q-1. 

[0057] In the definition of the decoder, let us assume that there is 
no interference in the output of the decoder. There are two reasons for this: 
firstly, this ensures that if the censoring has been completely successful, the 
5 remaining samples can be processed optimally. Secondly, if the interference 
signal is weak, the loss that formula (5) has caused in the signal, compared to 
the conventional decoder processing the original received signal, is minimised. 

[0058] The soft decision on the i:th Map decoder is obtained as a 
logarithm of the APT (a posteriori likelihood) ratio of each information bit dk = 1 
10 compared to APT = 0. Thus the soft decision has the form 

Zn^In.-^ i(Yh"(k).m',m) a,_i(rTT)^,(m)) 
' ' Z^Z^.ro(Yh''(k),m',m)ak-i(nT)ySk(m))" 

[0059] Forward and backward recursions of the Map decoder can 
be expressed as follows: 

I 

ak(m) = consta ZZ^iW (k),m'.m)«k-i(m') 

m' I— 0 
1 

15 Pk(m) = constp EEr.K(kXm",m);Sk_i(m') 

where consta ja constp are constants caused by normalisation. Transition 
probability yi() , where i = 0,1 , is, for its part, is obtained in the following 
manner. 

20 Yi{Yh (k),m,m') = Prob(Yh^(k) | dk=i, Sk=m. Sk-i=m')Prob(dk=i, Sk=m, | Sk.i=m' 
)■ 

Following the earlier mentioned publication of M. Lops and A. Tulino and its 
markings, the above mentioned formula for both type I and type II censoring 
25 can be expressed in the following form: 

Q 

Yi(Yh^(k),m,m") = Prob (dk=i) Prob (Sk+i=m' | Sk=m ) Z g(|vi|....,|vN|,q=q^"^) x 

^ ||(vP(k + (h-1)Re+j)-GP(k + (h-1)Re)UAej'*sf 
Qexpi! — L (9) 

j=0 -^"^o 
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where q^"^ is a variable having Q possible values representing the n:th 
censoring event where 0 = 2^^ for type I censoring and Q = N! for type II 
censoring. It is to be noted that the term g(|vi|,...,|vN|,q=q^"^) does not depend 
on the code bit b(l) and is clearly a constant term, which is removed in the 
5 reliability calculation. The first term of formula (9), a priori information on the 
transmitted information bits, is a function of the external Lext information. A 
reference is herein made to the publication P. Robertson, Illuminating the 
Structure of Code and Decoder of Parallel Concatenated Recursive 
Systematic (Turbo) Codes, Proc. Globecomm '94, December 1996. As 

10 appears from the presented formula (9), in connection with the metric 

calculation, filtering adapted to the spreading code, i.e. DS decomposition 
decoding, is also performed. This is indicated by the term s in formula (9). In 
order to operate the presented LogMap algorithm requires the information on 
the signal-to-noise ratio. If the estimate of the signal-to-noise ratio is not 

1 5 available, a so-called SOVA algorithm known to a person skilled in the art can 
be applied. The publication W. Feng, B. Vucetic, A List Bidirectional Soft 
Output Decoder of Turbo Codes, International Symposium on Turbo Codes, 
Brest, France, 1997 is given as an example of this. The branch metric for both 
censoring types is as follows: 

R-1 

20 A ^ (dk) = -2 X 9^ { Ae'^s^ U^^ v (k+ (h-1 )Re +j) } 

where ()^ refers to conjugate transpose. 

[0060] Let us examine an example of an embodiment of the 
invention by means of the flow chart of Figure 5a. A signal is received in step 
500. In step 502, a default value, for example 0, is determined for the 

25 estimates of the desired signal for the first round of the interference 

cancellation iteration. In step 504, the estimate is subtracted from the received 
signal. In step 506, an interference subspace is defined and its significance in 
the received signal is reduced. In step 508, turbo decoding is performed and 
the estimates for the desired signal are determined. These estimates are 

30 passed back to the interference cancellation 510. 

[0061] Flow chart 5b shows another example. The chart is 
otherwise similar to the previous one, but after the turbo decoding, a decision 
is made in step 512, whether a new iteration is performed for the decoding. If 
not, the estimates 510 are passed to the interference cancellation, if yes, the 

35 estimates 514 are passed to the input of the decoder again. 
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[0062] Although the invention has been described above with 
reference to the examples according to the attached drawings, it is obvious 
that the invention is not restricted thereto, but may be modified in a variety of 
ways within the scope of the inventive idea disclosed in the attached claims. 
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